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We study several antiferromagnetic formulations of the 0(3) spin model in three dimensions by means of Monte 
Carlo simulations. We discuss about the vacua properties and analyze the phase transitions. Using Finite Size 
Scaling analysis we conclude that all phase transitions found are of first order. 
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1. Introduction 

The study of antiferromagnetic models (AFM) 
may give us some insights on the formulation of 
non trivial relativistic quantum field theories, as 
they present a very rich phase space and, pre- 
sumably, new universality classes They are 
also interesting as related with models describing 
high temperature superconductors |Q ; the two di- 
mensional quantum Heisenberg antiferromagnet 
in the low temperature region can be described 
by a non linear a model in three dimensions Q . 

We study in this contribution several formu- 
lations of three dimensional antiferromagnetic 
sigma models as a first step in these directions. 

It is clear that we must go beyond the naive 
formulation of an 0(3) a model in a cubic lattice 
with only nearest neighbors coupling 
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where <Ti is a three components normalized real 
vector, as it can be easily seen that the antiferro- 
magnetic system (/3 < 0) is trivially related with 
the ferromagnetic one obtained by the staggered 
transformation 



(2) 



We will consider three models that break the 
symmetry (0) making the vacuum frustrated: 



1. A two parameter model with nearest (nn) 
and next to nearest (nnn) neighbor couplings 



S = (3l'^CTi ■ (Tj 
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Under the transformation (g) the second sum in 
(||) does not change so it is not possible to map 
the negative [32 values onto positive ones. 

2. A face centered cubic lattice (FCC) whose 
geometry explicitely breaks the symmetry (^. 
Other cubic lattices like the BCC (interior cen- 
tered cubic) and the tetrahedrical (diamond) fail 
to do so. Notice also that by setting /3i = in 
(||) two FCC sublattices are decoupled. 

3. A Fully Frustated model constructed by 
defining the following set of couplings 



(4) 



where (3x.y,z;^i is the coupling of the link pointing 
in the fi direction from the x, y, z lattice site, the 
values /i = 0, 1, 2 correspond to the x, y, z direc- 
tions respectively. This model presents a Z2 local 
gauge symmetry: it is invariant under a change of 
the sign of a particular spin and a simultaneous 
change of the sign of the couplings at the links 
starting from the same site. 
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2. The simulation 

We have used mainly the Metropolis algo- 
rithm for the updating with several overrelax- 
ation steps. For the largest lattice size (L = 64) 
the number of Monte Carlo sweeps performed af- 
ter thermalization, has been of the order of 10^, 
while for the smaller lattices that number has 
been even greater. We have checked in all cases 
that the autocorrelation time was much smaller 
than the total Monte Carlo time used for mea- 
sures. 

We have found the Wolff's single cluster algo- 
rithm Q to be very inefficient near to the anti- 
ferromagnetic transition, because the size of the 
biggest cluster ussually represents a very large 
fraction of the total lattice volume. For this rea- 
son, we have used it only to study the ferromag- 
netic transition. 

As observables we have measured the energies 
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with V being the lattice volume. In the case of 
the FCC lattice we only measure while in the 
Fully Frustrated model (see |^ for details) the 
sign of the coupling has to be properly taken into 
account. 

For antiferromagnetic phases the standard def- 
inition of the magnetization as M = y cTj is 
not an order parameter. For the first two models 
we have instead considered a staggered magneti- 
zation defined as 



(7) 



where a = x,y, z. For the third model we have 
constructed the following set of vectors 
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(even) 
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with i, j,k = 0,1. It can be checked that the mean 
values of the previous quantities are independent 
of i,j,k. We shall refer to it as the period two 
magnetization. 



In practice we measure the magnetizations 
squared from which we compute the Binder cu- 
mulant and the susceptibility defined respectively 



Ul = {M^)/3{M^) 
X = V({M')-{\M\[ 



(9) 



Another interesting quantity is the correlation 
length. To avoid problems with fluctuations and 
asymmetric lattices we use the second momentum 
definition considered in ref. Q , valid only for the 
disordered phase 



go/ 91 - 1 

4sin2(7r/L) 
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where go and gi are the Fourier transforms of 
the propagator G{ri — rj) — {cTi ■ a-j) at zero and 
minimal nonzero momentum respectively. 

The derivatives of the energies and magnetiza- 
tions with respect to the couplings can be com- 
puted as connected correlations. For instance, the 
specific heat matrix can be expresed as 

ri E 

= ^ = 3VmE,) - {E,){E,)) (11) 



We have used a Finite Size Scaling Analysis to 
compute the critical point and the critical expo- 
nents associated with the phase transitions. Mea- 
suring for instance the maxima of the specific heat 
and the susceptibility, using the spectral density 
method, we obtain 

C^L"'/" , x-L''^''- (12) 

The critical temperature can be obtained from 
the scaling behavior of the crossing point /^Li.La 

1 1 1 



Pc \og{L,/L2)- ^^^^ 

The critical exponent /3 may be computed from 
the magnetization 



/3c 



(14) 



For first order phase transitions however the scal- 
ing behavior presents fictitious critical exponents 
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3. Results 

We have analyzed the phase diagram of the 
two parameter model (see figure|l|). We have 
found three phases: ferromagnetic and antifer- 
romagnetic separated by a disordered (param- 
agnetic) phase. The order parameter for the 
paramagnetic-ferromagnetic (P-F) transition line 
is the magnetization, while for the paramagnetic- 
antiferromagnetic (P-A) we use the staggered 
magnetization. Along the P-F line we find, for 
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Figure 1. Phase diagram for the two parame- 
ter model. The solid line corresponds to the P-F 
transition and the dashed line to the P-A tran- 
sition. The point where /32 = is the standard 
cr model critical point. The points where /3i = 
correspond to the critical ferromagnetic and an- 
tiferromagnetic critical points of the FCC model. 



the following points in the (/3i,/32) plane 

A = (/3i ^ 2, /32 = -1.25111(13)) 

B = (/3i= 0.85763(8), /32 = -2/3i) (16) 

C = (/3i = 0.5, /32 = -2.3899(12)) 

We have found very difficult to reach the 
asymptotic behavior even for large lattice sizes. 
We should emphasize that it is even harder to at- 
tain at the point B. Assuming the results of the 
L = 64 as asymptotic we found that the growth 
of the specific heat is compatible with an a/i^ = 3 
while lower exponents can be readily discarded. 

From an analysis of the magnetic susceptibility 
we exclude a second order transition. The energy 
histograms (see figure^) confirm this asumption, 
for large enough L they present a clear two peak 
structure with an stable inter-peak distance. We 
point out that even for the smaller lattice sizes 
where the two peaks cannot be resolved the spe- 
cific heat and the susceptibility give strong indi- 
cations of the transition first order character. 
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(32 — the standard a model critical point, for 
/?i = the original lattice is decoupled into two 
sublattices and it is equivalent to the ferromag- 
netic FCC model with /? = 0.619(5). We have 
measured the critical exponents on those points 
along the P-F line, checking that they agree well 
with known values of the standard a model. The 
exponent v has been obtained with a 2% of accu- 
racy. 

Along the P-A line wc have measured along 
fixed directions the specific heat, susceptibility, 
staggered magnetization and correlation length at 



Figure 2. Energy histogram for L = 32 (solid), 
L = 48 (dashes) and L = 64 (dot-dashes) at B. 



To estimate the correlation length according to 
(p^), we have to simulate in the disordered phase 
and then extrapolate the results in order to obtain 
the values at the critical point. We have found 
the following values for the correlation lengths de- 
fined in (0) 

_ 7 _ ~ 7 , (17) 
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with statistical errors of the order of 10%. These 
values explain a posteriori the difficulties found in 
reaching, especially for point B, the asymptotic 
region. 

The second model considered is the FCC lat- 
tice with antiferromagnetic coupling. Classically 
the ground state presents a 0(3)^ degeneracy 
group. However, when thermal fluctuations are 
taken into account a coUinear ground state is se- 
lected (see for details). This is an example 
on Villain's order from disorder In the ab- 
sence of another interaction that fixes the global 
direction, there remains a Z2 degeneracy. 

We have performed a numerical simulation in 
the L = 24 lattice in the low temperature phase 
(/3 = —5 < (3c) that confirms the coUinear ground 
state structure. Each of the 2^ ground states is 
very stable under Monte Carlo evolution with a 
local update algorithm. Extrapolating to infinite 
volume the transition point, in the assumption of 
a first order behavior, we obtain 



(18) 



The energy histogram of the L = 32 lattice 
presents a clear double peak that establishes the 
first order nature of the transition, although for 
that lattice size the asymptotic behavior has not 
been reached. 

Finally let us present our results for the Fully 
Frustrated model. As the hamiltonian of this 
model is invariant under the transformation (|^) 
we will consider only the (3 >{) case. We obtain a 
phase transition, between a disordered phase for 
small (3 and an ordered one with a complicated 
structure, at 
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2.26331(13) 



(19) 



The ground state in the ordered phase is highly 
degenerate. We have checked for L = 2, by 
means of numerical and analytical methods, that 
the equilibrium configurations, in the {61,62,63} 
space where = cr(ro) • c(T'i), lay inside the 
hexagon perpendicular to the (1,1,1) vector with 
vertices at (1, -^) and permutations. For 

L > 2 the equilibrium configurations concentrate 
around the six corners of the hexagon in a region 
whose size decreases with increasing lattice size. 



The finite size scaling analysis of the specific 
heat shows that this is a weak first order phase 
transition with a hard to reach thermodynamic 
limit. From the L = 8 to L — 24 values we 
find good agreement with a/v = 1 while from 
the L = 24 to L = 64 we find an increasing value 
that tends to a/v = 2> as expected for a weak 
first order phase transition. We only observe an 
incipient two peak structure in the energy his- 
togram for L — 64. For the susceptibility we find 
a behavior compatible with a first order charac- 
ter. From the bigger lattices we quote ^ /v ~ 3. 
We therefore conclude that the order parameter 
for the transition, i.e. the period two magnetiza- 
tion becomes a discontinuous function of (3 in the 
thermodynamic limit. 

4. Conclusions 

We have explored three models with internal 
0(3) symmetry that are antiferromagnetic or de- 
velop frustration. We conclude that all the anti- 
ferromagnetic transitions found are of first order. 

We have checked that the ferromagnetic line 
belongs to what seems to be the only universality 
class for 0(3) models in three dimensions. 
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